Negative Power Spectra in Quantum Field Theory 
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We consider the spatial power spectra associated with fluctuations of quadratic operators in field 
theory, such as quantum stress tensor components. We show that the power spectrum can be 
negative, in contrast to most fluctuation phenomena where the Wiener-Khinchine theorem requires 
a positive power spectrum. We show why the usual argument for posit ivity fails in this case, and 
discuss the physical interpretation of negative power spectra. Possible applications to cosmology 
are discussed. 
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The well-known Wiener-Khinchine [TJ [2] theorem 
states that the Fourier transform of a correlation func- 
tion is a power spectrum. A corollary of this theorem is 
that the power spectrum can normally be written as the 
expectation value of a squared quantity, and hence must 
be positive. The theorem is most commonly formulated 
for temporal Fourier transforms. (For a review, see for 
example, Ref. [3 .) However, it may also be formulated 
for spatial Fourier transforms, which will be our focus. 
The purpose of this letter is to note that there is a loop- 
hole in the proof of the corollary, which allows for the 
possibility of negative power spectra. Furthermore, we 
show that such negative spectra can arise in the case of 
the fluctuations of quadratic operators in quantum field 
theory, such as the energy density operator. 

First, let us recall the more familiar case where the 
power spectrum must be non-negative. Let x) be a 
fluctuating quantity, e.g., a Hermitian quantum opera- 
tor in fiat spacetime, and let the associated correlation 
function be 

C(t-t',x-x') = (F(t,x)F(t',x')). (1) 
Define k) to be a spatial Fourier transform 

%k) = ^a Jd 3 xe ik -F(t,x), (2) 

and let its correlation function be 

C(t-t',k,k') = (F(t,k)F(t'X)). (3) 

If we use insert Eq. Q into Eq. (J3|, and change integra- 
tion variables to u = x — x' and v = x + x r , we find 

C(t-1?,ky) = - 7 ^5(k-k!) [ dW^C^-^u). 
4(27r) d J 

(4) 



Hence 

C(0,k,k') = P(fc) £(k - k') , (5) 
where the power spectrum P(k) is defined by 



p(k) 



(2tt) 3 



tfW ku C(0,u). 



(6) 



If C(0,k, k') exists, then it is non-negative, as it is 
the expectation value of F 2 (£,k), and as a consequence 
P(k) > 0. This is a form of the Wiener-Khinchine theo- 
rem for spatial Fourier transforms. 

A possible loophole can arise if C(0,k, k') is not well- 
defined. A simple example can illustrate this possibility. 
Consider the operator 



F(*,x) =:y> 2 :, 



(7) 



where cp is a massless free scalar field in four-dimensional 
Minkowski spacetime. The vacuum correlation function 
in coordinate space may be found from Wick's theorem 
to be 



C(T,r) 
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8tt 4 (r 2 - r 2 ) 2 ' 
where r = t — t'. The spatial Fourier transform is 



(8) 



CM) = ^/<iW-C M = -gg. (9) 

This result may be verified either by direct evaluation of 
the above integral using contour integration methods, or 
by checking that 



C(r,r) = J d 3 fce- ik,x C(r,fc) 
1 f°° 

— A — lim / dkk sin(/cr) sm(kr)e~ ak .(10) 
167r 4 rr a^o J K J 
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The associated power spectrum is 
P(k) = (7(0, k) = - 



k 
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(11) 



which is negative. Clearly, the positivity argument given 
above does not hold in this case, and the reason must be 
that C(0,k, k') does not exist as a well-defined quantity. 

The explicit form for F(x, t) in terms of creation and 
annihilation operators is 



iKq-qO-x-C^-u/)*] t 



f(*, X ) = y — l -^U 

^ 2VvW 1 
q,q 

+ e i ^ q+q/ > x -^ +u;/ ^ a q a q/ + #.c.} . (12) 

Here we use box normalization in a volume V, H.c. de- 
notes the Hermitian conjugate, = |q|, and Co>' = |q'|. 
The Fourier transformed operator is 



F(t,k) = 



(Wi-£j)t t 




x q+k u q 



(13) 



where cji = |q + k|. 

The correlation function for F(k, t) can be shown to 

be 

where k = |k|. If we take the limit r — > 0, this quantity 
is formally positive, but divergent and hence undefined. 
This is the positive quantity that would appear in the 
Wiener-Khinchine theorem and implied a positive power 
spectrum had it been well-defined. If we first take the 
large V limit and perform the integration on q, then the 
result is 



C(r,k,k') 



where we have used V ^k,k / /(^ 7r ) 3 S(k — k'). Now we 
can take the r — » limit, and obtain the negative power 



spectrum in Eq. (11). 



A more physically interesting example arises in the case 
of the fluctuations of quantum stress tensors, such as the 
energy density of the electromagnetic field. Let p(t,x) 
be the normal-ordered energy density operator, and its 
vacuum correlation function be 



C (r,r) = (p(i,x)p(i / ,x')} = 



(T 2 +3r 2 )(r 2 + 3r 2 ) 



7r 4 (r 2 - r 2 ) 6 

The explicit form of its spatial Fourier transform is 

fc 4 sin(A;T) 



C {r,k) 
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(16) 
(17) 



P(k) 
\Po\ 




FIG. 1: The vacuum and thermal power spectra, Po and Pr, 
for the electromagnetic field are plotted as functions of tem- 
perature. Here we use units in which Boltzmann's constant 
is unity, ks — 1. The total power spectrum Po + Pt, is also 
plotted, and seen to be positive for T > 1.04/c. 



Again, this result may be confirmed either by contour in- 
tegration, or by checking that the inverse Fourier trans- 
form of Co(t, k) is indeed Co(r, r), as illustrated in 
Eq. (10). The power spectrum is 



P (k) = C (0,k) 



k 5 
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(18) 



which is again negative. 

This leads us to the question of the physical interpre- 
tation of negative power spectra. One possibility is that 
the vacuum fluctuation contribution is a sub-dominant 
contribution which reduces the magnitude of a net posi- 
tive power spectrum. This possibility can be illustrated 
by energy density fluctuations at finite temperature. The 
electromagnetic energy density correlation function at fi- 



nite temperature can be obtained from Eq. (16) as an im- 
age sum over imaginary time by replacing r by r + in/3, 
and summing over all integers n. Here j3 = 1/(/cbT), 
where ks is Boltzmann's constant, and T is the tempera- 
ture. The temperature dependent part of the correlation 
function can be written as 



C T (0,r) 



71=1 

2 ^ (3r 2 -n 2 |8 2 )(r 2 -3n 2 ^ 2 ) 



7T 



n=l 



(r 2 +n 2 /3 2 ) 



when At = 0. The spatial Fourier transform of this ex- 
pression is Pr(k), the temperature dependent part of the 
power spectrum, which is found to be 



Pr(k) 



k 4 
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ln(l 



(20) 



Note that Pr(k) > for all fc, as may be seen in Fig. [I] 
We also see from this figure that the total power, Pq+Pt, 
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is positive for T > 1.04k /kp- Thus at high tempera- 
ture, the negative contribution from the vacuum fluctu- 
ations reduces the magnitude of the net power. At lower 
temperatures, the vacuum term dominates, and the total 
power is negative. 

Consequently, we also need a physical interpretation of 
the case of net negative power. Note that negative power 
spectra are always associated with coordinate space cor- 
relation functions which are singular at coincident points, 
and hence cannot represent the expectation value of a 
meaningful squared quantity. However, the correlation 
function at distinct points is meaningful, and can have 
either sign. If C(t — £',x — x') > 0, the fluctuations at 
(t,x) are correlated with those at (t^x'), and conversely 
if C(t — t' , x — x') < 0, they are anticorrelated. Changing 
the sign of the power spectrum, P(k) — >• — P(fc), changes 
the sign of C(t — t\ x — x'), and hence interchanges cor- 
relations and anticorrelations. So long as (7(0, 0) is unde- 
fined, both situations are logical possibilities. Consider 
the contribution of a finite bandwidth to a spatial corre- 
lation function and define 



Ak 



(0,r)= / 

Jk, 



d d ke- lkx P(k). (21) 



ko<k<ki 



This quantity would describe the spatial correlations in a 
situation where fluctuations with k < ko or k > k\ have 
essentially been filtered out, It will typically be a quasi- 
oscillatory function in space. Changing the sign of the 
power spectrum in this interval interchanges the minima 
and maxima of Ca&(0, r), interchanging correlations and 
anticorrelations. 

In general, it is not C(t — t' ,x — x') itself, but rather 
integrals of the correlation function over finite spacetime 
regions which are observable. Let S\{x) and S2(x') be 
test functions which describes the effect of some measur- 
ing apparatus. The correlation function for the outcomes 
of the measurements described by Si and S 2 is 



d^xS^x) / d i x'S 2 (x') C(£-£',x-x'). (22) 



K 



Even though the function C is singular at coincident 
points, it is well-defined as a distribution, so K is finite. 
Consider the case of the vacuum energy density of the 
electromagnetic field, where C is given by Eq. (16). This 
correlation function may be expressed as a total deriva- 
tive |4j: 



0)(r,r) = - 



1 
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v 2 nv u n' hS[(x-x'y/e 2 ], (23) 



where £ is an arbitrary length. Here V 2 and □ denote 
the Laplacian and d'Alembertian operators, respectively, 
in x, and V' 2 and the corresponding operators in x' . 
The value of Co(t, r) is unchanged if £ changes. We next 
integrate by parts in the expression for K and assume 
that the surface terms vanish, which will be the case if S\ 



and S 2 vanish sufficiently rapidly at infinity. The result 
is 

K = J AV 2 Dft(x) 

x J d 4 x'V 2 n' S 2 (x') \n 2 [(x-x , ) 2 /£ 2 }. (24) 

This expression contains only an integrable, logarithmic 
singularity and is hence finite. [An explicit example for 
K is illustrated in Fig. 5 in Ref. [4].] The key point is 
that the singular nature of the coordinate space correla- 
tion function, which always accompanies negative power 
spectra, does not prevent quantities such as K from being 
well-defined. 

Radiation pressure fluctuations on a mirror can be 
computed as integrals of a quantum stress tensor cor- 
relation function [5 , yielding a result which may also be 
derived by an alternative approach based on photon num- 
ber fluctuations [6] . This is an illustration of how prop- 
erly defined integrals of singular correlation functions are 
physically meaningful. 

Quantum stress tensor fluctuations can potentially 
contribute to the primordial density fluctuation spectra 
in inflationary models [7J [8] . Here the observable quan- 
tities involve time integrals of the energy density cor- 
relation function. Even though the integrands are sin- 
gular, the integrals are finite. In a model with a sin- 
gle scalar inflaton field, the contribution of electromag- 
netic energy density fluctuations to the density fluctua- 
tion power spectrum can be written as [See Eq. (88) of 
Ref. El.l 



Ps P (k) 



£ P HS 2 
30 tt2 



4ttH 
bk 



(25) 



where £p is the Planck length, H is the Hubble param- 
eter during inflation, and S is the scale factor change 
during inflation. Note that the first term is negative. In 
the models discussed in Ref. [8] , this term is dropped, as 
it corresponds to a delta-function term in the coordinate 
space correlation function, and is hence not observable 
in measurements made in disjoint regions. However, if 
there is any process which has the effect of filtering this 
negative power spectrum, it would no longer be a delta- 
function and hence could become observable. Classical, 
non-linear evolution of the density perturbations is one 
possible filtering mechanism. Thus the appearance of this 
term illustrates the possibility of negative power spectra 
in cosmology. The appearance of a negative power spec- 
trum of tensor perturbations in inflationary cosmology 
will be discussed in a separate publication [9]. 

We have focused on fluctuations in space, and power 
spectra defined by spatial Fourier transforms. As noted 
earlier, the Wiener-Khinchine theorem can also be for- 
mulated in terms of temporal Fourier transforms and an 
analogous power spectrum P(uS) can be defined. If the 
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correlation function in time is finite in the coincidence 
limit, then P(u) > 0. In the case of examples studied 
earlier, the coincidence limit does not exist, but nonethe- 
less one finds a positive power spectrum. For example, 
for the case of the vacuum fluctuations of the electromag- 
netic energy density, the power spectrum for temporal 
fluctuations is 



P(w) 



dre 



C o (r,0) 
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> 0. 



(26) 



It is not known whether there are examples of negative 
temporal power spectra. 

In this letter, we have illustrated how the fluctuations 
of quadratic operators in field theory can produce neg- 
ative power spectra for spatial fluctuations. A negative 
contribution can have the effect of decreasing a net pos- 
itive spectrum. However, it is also possible for the net 
spectrum to be negative, as in the case of vacuum en- 
ergy density fluctuations of the quantized electromag- 
netic field. Negative power spectra have the opposite 
correlation-anticorrelation behavior as does a positive 
spectrum with the same functional form. 
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